Nonequilibrium thermodynamics of a squeezed harmonic oscillator 
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We consider the thermodynamic properties of the squeezed vacuum state of a frequency- 
modulated quantum harmonic oscillator. We analytically relate the squeezing parameter to the 
irreversible work and the degree of nonadiabaticity of the frequency transformation. We furthermore 
determine the optimal modulation that leads to maximal squeezing, and discuss its implementation 
as well as the detection of squeezing in single cold ion traps. 
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The time-dependent quantum harmonic oscillator 
serves as a model system for a variety of physical prob- 
lems. Notable examples include the generation of non- 
classical states [ij and the dynamics of cold ions in Paul 
traps [H. More recently, it has played a major role in the 
study of cosmological particle creation and nonequi- 
librium quantum thermodynamics An experimen- 
tal investigation of a harmonic atom trap with frequency 
jumps has moreover been reported in Ref. The im- 
portance of the time-dependent quantum oscillator stems 
from the fact that it is an exactly solvable system. The 
propagator and the transitions probabilities of an oscilla- 
tor with time-dependent frequency and time-dependent 
linear driving have been first derived by Husimi using 
a Gaussian wave function ansatz Q. A more general 
method based on invariant operators has been later de- 
veloped by Lewis and Riesenfeld 0. 

A remarkable property is that the state of the oscil- 
lator is squeezed when its angular frequency is changed 
in a nonadiabatic way [1, The degree of squeezing 
depends on the specific frequency modulation considered 
and various protocols have been discussed [Tol. [ill. [l^. ITst . 
Efficient schemes for the production of squeezed states 
are nowadays essential for high precision measurements, 
such as in interferometric gravity wave detectors [l^ . In 
addition, squeezing appears as a crucial resource for en- 
tanglement creation in quantum information applications 
with continuous variables [Tsl . [T^l . The determination of 
an optimal squeezing protocol is thus of high importance. 

The quantum harmonic oscillator with arbitrary fre- 
quency modulation has recently been studied from a ther- 
modynamical point of view in Ref. ^17] . Using the gen- 
eral formula for the probability density of quantum work 
introduced in Ref. [l^, the statistics of the total work 
performed on an isolated, but initially thermal, oscillator 
during adiabatic and nonadiabatic variations of its angu- 
lar frequency has been determined. In particular, an ex- 
pression for the irreversible work, also often referred to as 
dissipative work [lo'l, has been obtained. The irreversible 
work is defined as the difference between the total work 
{W) and the free energy difference AF (the reversible 
work): {Wirr) = (W") - AF. The total work is here 
given by the difference between final and initial energies 
of the oscillator, (VF) = {H{t)) - {H{0)), where H{t) is 
the time-dependent Hamilton operator the quantum os- 



cillator. On the other hand, the free energy difference can 
be written in terms of the partition function Z{t) in the 
usual manner as AF ^ F{t)-F{0) = -kT In Z{t)/Z{0), 
with T the initial temperature of the oscillator. The ir- 
reversible work corresponds to the energy absorbed by 
the system during a nonadiabatic transformation and is 
therefore zero for an adiabatic, reversible, change of the 
frequency. It also represents the energy that would be 
dissipated into a heat bath at temperature T, were the 
system weakly coupled to one [lOl . The irreversible work 
has been related to a parameter Q* (see Eq. (ITU)) be- 
low), originally introduced by Husimi. The latter can 
be regarded as a measure of the degree of nonadiabatic- 
ity [17]: Q* is unity for an adiabatic transformation of 
the oscillator's frequency and increases monotonically the 
less adiabatic a transformation is. 

In this paper, we express the amount of squeezing gen- 
erated by an arbitrary frequency change in terms of the 
nonadiabaticity parameter Q* and the irreversible work 
(Wirr) for an oscillator initially in the ground state. By 
considering the thermodynamics of vacuum squeezing 
and, in particular, by quantifying the degree of squeezing 
with the help of the nonadiabaticity parameter, we are 
able to extend initial studies of squeezing production in 
frequency-modulated harmonic oscillators. We further 
use optimal control theory [2l| to determine the modu- 
lation that leads to maximal squeezing and compare the 
results with those obtained in Ref. [l^l . We finally discuss 
the experimental implementation of the optimal modula- 
tion as well as the determination of the thermodynamic 
and squeezing properties of the quantum oscillator using 
single ions in linear Paul traps. 

Thermodynamics of squeezing. A quantum harmonic 
oscillator with time-dependent frequency Lu{t) and unit 
mass is described by the Hamilton operator. 



(1) 



We consider a protocol where the frequency is changed 
from an initial value uj(0) — loq to a final value i-u{t) — ui 
during time r. The solution of the Heisenberg equations 
of motion for the position and momentum operators q{t) 
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and p{t) can then be written in the general form, 

q{t) = q{0)Y{t)+p{0)X{t) , (2) 
Pit) = m , (3) 

where the two functions X{t) and Y{t) are the solutions 
of the equation of motion of the corresponding classi- 
cal oscillator, X + LLj'^{t)X ~ 0, with initial conditions, 
X{0) = 0, X{0) = 1 and Y{0) = 1, Y{0) = 0. The 
latter ensure that the canonical commutation relation, 
[x,p] = i, is satisfied (we set ^ = 1 throughout). The po- 
sition and momentum variances for an oscillator initially 
in the ground state follow accordingly as. 



(qp) = 



Y^ 



ujoX^ 
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(4) 
(5) 
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As shown by Husimi, details about a specific protocol 
uj{t) is fully contained in a function Q*{t) defined as 0: 



Q*it) = 



1 



uj^{uj^{t)X^ + X^) + uj^{t)Y^ + Y' 



By using Eqs. ^ and ([51), we easily find that 
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where the parenthesis is recognized as twice the expecta- 
tion value of Hamiltonian ([T]). We thus obtain that the 
mean energy of the oscillator at any given time is simply 



{Hit)) 



Q*it) 



(9) 



The above equation provides some insight into the phys- 
ical meaning of the parameter For an adiabatic 
transformation, Q*it) = 1, and the mean energy of the 
oscillator is just given by the ground state energy. This 
corresponds to the known classical result that the ac- 
tion of the system, defined as the ratio of the energy 
and the angular frequency, is a time-independent con- 
stant. On the other hand, for a nonadiabatic change of 
the frequency, the value of > 1 indicates how far 
the final state of the oscillator is from its initial (equilib- 
rium) ground state. The latter statement can be made 
more precise by directly computing the irreversible work 
(Wirr) at the final time t. By evaluating the total work 
(W) and the free energy difference for Hamiltonian 
(P, one finds [ij, 



{W„r) 



coir) 



iQ*ir)-l) 



(10) 



The irreversible work is therefore zero for adiabatic trans- 
formations, as expected, and grows linearly with the 
nonadiabaticity parameter Q*it). 

Let us now establish a relationship between the irre- 
versible work and the degree of squeezing of the harmonic 
oscillator. At any given time, a squeezed oscillator state 
can be parameterized as (22l |. 
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sin^e + e^''cos^6 



sinh(2r) sin 9 cos 9 . 



(11) 
(12) 
(13) 



The time dependence of the squeezing parameter r{t) 
and the rotation angle 9{t) is controlled by the frequency 
modulation Lo{t). One should note that the mean values 
(x) and (p) remain here zero at all times. By inserting 
Eqs. pT|) - (fT^ into Eq. ([H]), we obtain the following rela- 
tion between Q* and the squeezing parameter r, 



Q*{t) = cosh2r(i) . 



(14) 



Equation p4|) is an important result that directly con- 
nects the degree of squeezing of the oscillator to the nona- 
diabaticity parameter, and thus to the frequency modula- 
tion Lu(t). It clearly shows that squeezing requires a nona- 
diabatic change of the frequency, Q*(i) > 1, and that 
large squeezing implies an exponential increase of Q*(i): 
Q*{t) - exp(2r(i))/2. By further combining Eqs. dll 
and (1141). we arrive at 



{W^rr) = ^^(t) siuh^ r(r) 



(15) 



(8) Two limiting cases follow from this equation: (Wirr) 



UJ r for small squeezing and (Wirr) ^ uJ e /4 for large 
squeezing. A high degree of squeezing thus requires an 
exponentially large amount of irreversible work. Equa- 
tions (fTO| . (fT4|) and (fT5|) reveal the intimate and simple 
relationship existing between the vacuum squeezing prop- 
erties of the time-dependent harmonic oscillator on the 
one hand and its nonequilibrium thermodynamic proper- 
ties on the other: the knowledge of the squeezing parame- 
ter allows the determination of both the degree of nonadi- 
abaticity of the frequency modulation and the amount of 
irreversible work produced. Conversely, the knowledge of 
the nonequilibrium thermodynamic state of the quantum 
oscillator gives direct access to its vacuum squeezing. 

Optimization of squeezing. An important question 
from a theoretical as well as practical point of view is 
the determination of a frequency protocol that leads to 
maximum squeezing for a prescribed maximal modula- 
tion amplitude. We use in the following optimal con- 
trol theory to answer this question for fixed values of the 
initial and final frequencies loq and toi of the harmonic 
oscillator. Optimal control theory (OCT) is a powerful 
mathematical optimization method based on the calculus 
of variations [21 1 . It allows to determine the function that 
minimizes a given cost functional in analogy to the famil- 
iar Euler-Lagrange equations of classical mechanics that 
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FIG. 1; (color online) Frequency modulation of the harmonic 
oscillator as a function of time for fixed initial and final fre- 
quencies Luo — 1 and uji = 1.3. The dashed (red) line is the 
result obtained by optimal control theory and the solid (blue) 
line corresponds to the result of Ref. [I4]. 



minimize the aetion of a system. The results of a numer- 
ical implementation of OCT using Pontryagin's principle 
for Wo = 1 and wi = 2 are summarized in Figs. 1, 2 and 3. 
The oscillator is taken to be initially in the ground state 
and the cost functional to minimize is chosen as 1/ {H{t)), 
since the average energy is a monotonous function of the 
squeezing (see Eqs. (0) and (O). Figure 1 shows the 
optimal frequency modulation as a function of time; we 
observe that it consists of a regular sequence of frequency 
jumps at which u}{t) abruptly switches from uj^ to uji and 
back. The corresponding squeezing parameter, depicted 
in Fig. 2, increases by discrete increments at each of the 
frequency jumps until the modulation is terminated. 

It is interesting to compare the above frequency pro- 
tocol obtained with the help of OCT to the one stud- 
ied analytically by Janszky and Adam 



13. Although 



the two protocols are not exactly identical, they lead to 
very similar squeezing results (see Figs. 1-3). Janszky 
and Adam considered a sequence of sudden frequency 
changes between loq and wi, separated by some delay 
times Ti . They found that squeezing was strongest when 
these delay times were given by a quarter of the oscilla- 
tion periods of the oscillator, that is, tq = 7r/2aJo before 
a jump from ujq to uii and ti = 7r/2ci;i before a jump 
from (jji to (jJo- The latter condition exactly corresponds 
to the time needed to exchange position and momentum 
axes in phase space. The total squeezing generated after 
n such cycles was shown to be 



o2r 



(16) 



The degree of squeezing achieved by such a protocol 
therefore increases exponentially with the number of cy- 
cles. It then follows from Eq. (fT4|) that the nonadiabatic- 
ity parameter also grows exponentially with n for large 
squeezing, Q*{t) ^ {wi/wq)^ /2. We have checked that 
the frequency modulation of Janszky and Adam is actu- 
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FIG. 2: (color online) Vacuum squeezing of the harmonic 
oscillator generated by the frequency modulation shown in 
Fig. 1. The dashed (red) line is the result obtained by op- 
timal control theory and the solid (blue) line corresponds to 
the resuh of Ref. [13. The dotted ( green) line shows the non- 
optimal squeezing produced by a periodic driving at frequency 
2ljo, <^(i) = li'o + (i^i — <^o) sin(2ajoi)/2. 



ally a stable solution of the optimal control algorithm. 

Implementation in ions traps. A squeezed vacuum 
state of a quantum harmonic oscillator was created and 
observed experimentally using a single ion confined in a 
Paul trap [23] . Paul traps are ultrastable rf traps that al- 
low to prepare, manipulate and measure quantum states 
with high precision j24] . The first observation of quantum 
jumps was performed in a single ion trap [25l |. The fre- 
quency of a trap is determined by the external electrode 
voltages and by the size of the trap. A voltage mod- 
ulation therefore directly leads to a modulation of the 
motional frequency. Since commercial electronic com- 
ponents can achieve fast voltage switching rates, as com- 
pared to the trap frequency, tailored frequency variations 
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FIG. 3; (color online) Nonadiabaticity parameter Q*{t), 
Eq. p4p . generated by the frequency modulation shown in 
Fig. 1. The dashed (red) line is the result obtained by op- 
timal control theory and the solid (blue) line corresponds to 
the result of Ref. [12j . 
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can be implemented. 

In the experiment (23j . the squeezed vacuum state of 
a harmonicaUy confined ^Be"*" ion was generated by cool- 
ing the ion to its motional ground state using sideband 
cooling, and by irradiating it with two Raman beams 
with a frequency difference 2u!q. The latter induces a 
parametric driving at frequency 2u!q, which squeezes the 
ground state of the ion. The squeezed state was detected 
by laser-coupling motional and electronic levels of the ion 
and observing the fluorescence signal of the ground state. 
The probability Pg{t) that the ion remains in the elec- 
tronic ground state after a given coupling time t depends 
on the level population of the motional degree of freedom. 
The probability distribution of the motional Fock state 
Pn is a known function of the squeezing parameter and 
is given by the Fourier transform of the time signal Pg (t) 
(mapped by sequential experimental runs with different 
final times). The squeezing parameter f3 — exp(2r) could 
then be deduced by fitting the function Pg{t), yielding a 
value (3 — 40. We can evaluate the corresponding nona- 
diabaticity parameter from Eq. to be Q*{t) = 20. 
The squeezing protocol used in the experiment is not op- 
timal and similar (or higher) squeezing values could be 
obtained in less time by employing the optimal squeezing 
modulation discussed in the previous section (see Fig.[2|). 



It is worth noticing that the squeezing parameter can 
also be determined from a measurement of the mean en- 
ergy of the oscillator, by combining Eqs. ([9]) and (|14|1 . 
The average energy of the oscillator at a given time can 
be evaluated from the measured level population P„ via 
the simple expression {H[t)) = ^„ fttj„(n + 1/2)P„. By 
proceeding this way, the squeezing parameter can hence 
be determined directly without doing any numerical fit. 

Conclusion. We have presented a relationship between 
the degree of nonadiabaticity, the irreversible work and 
the vacuum squeezing of frequency-modulated quantum 
harmonic oscillator. We have found that both the nona- 
diabaticity parameter Q* and the amount of irreversible 
work grow exponentially with large squeezing. We have 
moreover determined the optimal modulation that leads 
to maximal squeezing using optimal control theory and 
found that the result is very similar to the protocol in- 
vestigated analytically by Jansky and Adam 12]. We 
have in addition discussed the experimental implemen- 
tation in single ion traps and proposed a new method 
to measure the degree of squeezing and determine the 
nonadiabaticity parameter. 

This work was supported by the Emmy Noether Pro- 
gram of the DFG (Contract LU1382/1-1) and the cluster 
of excellence Nanosystems Initiative Munich (NIM) . 
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